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ABSTRACT 

This paper presents a qualitative evaluation of a communication satellite driven by a white noise 
sequence, and proposes a filtering solution to predict the dynamics of the satellite. The prediction of the 
position and velocity of the communication satellite though being perturbed by the effect of solar array, 
gravity waves and drag resistance, was achieved using a least square algorithm based filter to determine 
an unbiased estimate. The dynamics of the satellite assumed linear, were examined from three broad 
classifications viz a viz: Parametric estimation (Least Mean Square method), Noise in communication 
system and the Kalman Filter. The significant practical benefits of using the MATLAB program for the 
Kalman filter blockset in the simulations were established. 

KEYWORDS: Communication Satellite, White Noise, Least Mean Square, Unbiased estimates, 
Kalman filter 

INTRODUCTION 

Filtering is desirable in many situations in engineering and embedded systems. Radio communication signals are 
often corrupted or heavily correlated with noise. A good filtering algorithm can remove the noise from 
electromagnetic signals while retaining the useful information (Simon, 2001). For many space applications 
involving data communication, a large number of satellites have been launched into the earth orbit, for 
effectiveness in their operation, the attitude must be observable and the dynamics predictable. It is possible to 
obtain a time history of the attitude of a communication satellite by suitable instrumentation and telemetry 
(Ekejiuba, 1986). A near-earth satellite, for example, orbiting in the attitude range of 150 km to 450 km 
encounters small but non-negligible aerodynamic forces due essentially to the characterization of drag 
resistance, gravity waves in stably stratified atmosphere especially affects free and forced rotation of 
astronomical satellites (Sohn, 1995). The influence of major environmental forces on the attitude response of 
gravity gradient satellites using both analytical and numerical techniques presents a problem of major interest to 
the telecommunication engineer, This is a game theory approach; Nature tries to maximize the estimation error, 
the engineer tries to minimize the error (Simon, 2008). This work models the dynamics of a communication 
satellite driven by a white noise sequence, and proposes a finite-time signal filtering solution using the Kalman 
filter, in which state estimation was addressed. Theoretically, the Kalman filter is an estimator for what is called 
the linear-quadratic problem, which is the problem of estimating the instantaneous "state" of a linear dynamic 
system perturbed by white noise — by using measurements linearly related to the state but corrupted by white 
noise. The resulting estimator is statistically optimal with respect to any quadratic function of estimation error. 
Practically, the Kalman filter is one of the greatest discoveries in the history of statistical estimation theory and 
possibly the greatest discovery in the twentieth century. It has enabled humankind to do many things that could 
not have been done without it, and it has become as indispensable as silicon in the makeup of many electronic 
systems. Its most immediate applications have been for the control of complex dynamic systems such as 
continuous manufacturing processes (Grewal and Andrews, 2008). The Kalman filter provides an efficient 
computational means to estimate the state of a process, in a way that minimizes the mean of the squared error. It 
supports estimations of past, present, and even future states, and it can do so even when the precise nature of the 
modeled system is unknown. Our approach determines the best estimate: the measured values of the input and 
output of the system are required, where by the measured error covariance and the estimated error covariance 
must be approximately same or unbiased. 

Statement of the Problem 

[CGWIC 2007]The NigComSat-1 was the first African geosynchronous communication satellite , when it was 
launched into orbit, aboard a Chinese Long March 3B carrier rocket , from the Xichang Satellite Launch Centre 
in China. The satellite, which is the third Nigerian satellite to be placed into orbit , was launched into a 
geosynchronous transfer orbit and subsequently it was successfully inserted into a geosynchronous orbit , 
positioned at 42.5°E on 13 May 2007. It had a launch mass of 5,150 kg, and had an expected service life of 15 
year. The satellite, designed to operate in Africa, parts of the Middle East and southern Europe, is to serve 
mobile phone users, facilitate internet performance in remote regions and improve overall 
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telecommunications. The spacecraft was operated by NigComSat and the Nigerian Space Agency, NASRDA , it 
is been monitored and tracked by a ground station built in Nigeria's capital Abuja .On 10 November 2008 (0900 
GMT), the satellite was reportedly switched off for analysis and to avoid a possible collision with other 
satellites. According to Nigerian Communications Satellite Limited, it was put into "emergency mode operation 
in order to effect mitigation and repairs "( AFP . 2008). On November 11, 2008, NigComSat-1 failed in orbit after 
running out of power due to an anomaly in its solar array. The need to predict and achieve an unbiased estimate 
of the actual path and position of a communication satellite is of great concern, hence this study. 

Significance of Study 

This study will no doubt be of immense benefit in the tracking a communication satellite as it helps to predict 
the present, and future position of the satellite at any required time. Since received signals are heavily correlated 
with noise, filtering is desirable in other to control a dynamic system, for this application, it is not always 
possible or desirable to measure every variable that you want to control, Kalman filter is used for predicting the 
likely future courses of dynamic systems that people are not likely to control, it provides a means for inferring 
the missing information from indirect and noisy measurements also giving the best estimate of the signal output 
thereby improving Signal to Noise ratio and Efficiency which are the desires of the telecommunication 
Engineer. 

The State Estimation Problem 

The procedure employed in least squares estimation of a process is usually carried out for a sequence of 
different order process and noise models until the best and simplest possible model is obtained. The more 
general estimation problem can be formulated on the basis of maximum likelihood and Bayesian techniques 
(Athans, 1971; Kalman, 1960) using statistical information in terms of joint probability distribution functions. 
However, for the linear dynamic system with additive, zero-mean white Gaussian measurement noise defined in 
terms of mean values and variances, which will be appropriate for many practical problems, the least-squares 
solution formulated as a deterministic problem with appropriate weighting leads to the maximum likelihood 
estimate. (Oyediran, 2010) Computation of the optimal estimates should consequently rely on convergence of 
the iteration employed; this is accomplished through sequential filtration of the estimate. 

Here, we consider a time-series model of the form 

m m 

x{j) + Y d a i x{j -i) = Y< b i u O - f ) (!) 

i=\ i=0 

The partial fraction expansion required for obtaining a time solution to (1) is defined for distinct roots using the 
method of residues, as 

v=l k=\ \S-(T V ) 
where OC l ...OC k are the poles of the function Y(s) with 
multiplicities m l ...m k . 

The inversion integral was then written in the form 

y{t) = — f Y(s)e s, ds = V Y-^*-f*-y v , t > (3) 



so that the function Y(s) with contaminated noise becomes 



Y{s) = 



U(s) + white noise (4) 



where CC { , //. are the parameters of the system, and n is the order of the system 
Using the shift operator Z~ defined by Z x(j) = x( j — 1) , equation (1) becomes 

A(z ) A(z ) 
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where y(k) is the observed output signal, u(k) is the applied input signal, N is the number of samples and 
e{t ) is the noise sequence. The polynomial operators A, B, and C are defined as follows 






i=\ 



iyn 

Biz- 1 ) = £^z-'' 



i=\ 

NC 



C(z- 1 ) = l + Y,c' i z~ i (6) 



(7) 



If the continuous model in (4) is discretized we shall obtain a discrete time model of the form 

y(*) = E /' z _ r u(k) + Mk) 

where the parameter set a, b is related to the parameter set CC , f3 via the relation 

a i =-exp(-a i T) 

P 
and b t = — (l-expH^T)) (8) 

a. 

where T is the sampling interval in seconds. The infinite noise of the discrete observations due to aliasing that 
may result from the above discretization process is negligible. This is justified if u(k) and e(s) are 
independent of all k and s. This is a reasonable assumption as long as the identification is performed for data 
acquired from experiments, where u(k) is a priori known sequence. In some practical situations, this 

assumption is often violated when operating records are used because in such a case the input may depend on 
the output through feedback. 

The canonical model in (5) was made equivalent to the discrete-time model of (7) by satisfying the conditions 
(i) C\ = a\ (9) 



(ii) 



A( Z - 1 ) it\ + a iZ - v 



The statistical method of maximum likelihood was employed to optimize the probability of obtaining the 
expected result. Consequently, a loss function V (0) was defined as follows 

V(6) = ^£ 2 (k) (10) 

which will be minimized with respect to the system parameter set = \fC , b\ . The residues were defined by 

£(k) = y(k)-?^-±u(k) (11) 

Mz ) 

so that the values of the parameter set a' and b' that make V(0) in (10) minimum will be the estimates of the 
parameters of the system. The approach considered is one of finding the coefficients of the prediction model 

C{z ) c\z ) 



y 



so that the mean square prediction error 

f k 



V{e) = Y j y{k)-y -A- =2> 2 (*) (13) 

k=\ V^ V k=l 
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is as small as possible. By doing so the assumption of Gaussian distribution of the noise sequence £\k ) may be 
relaxed. In the model given in (7), the residues were obtained as 

( n b z -l \ 

e{k)=y{k)- £— ' T u{t) (14) 

For system order greater than n, £\k) as given in (14) becomes computationally difficult and highly nonlinear, 
whereas in the model of (7), the residues given by (11) may be computed quite easily for any given system 
order. Since the two models are equivalent via the transformations defined in (9), we are free to use any of them. 
Thus the form of residues given in (11) is recommended and used in this study for parameter estimation, and 
application of a filter for signal noise attenuation becomes evident. 



Measurement Noise Filtration 

The relationship between Information, Bandwidth and Noise 
Rate 

R= 2Blog 2 (-ra) b/s 

Signal to Noise Ratio 

SNR=10log w (S/N) dB 



(15) 



(16) 



(Jianfeng, 2007) There is a theoretical maximum to the rate at which information passes error free over a 
channel. This maximum is called the channel capacity, C. The famous Hartley- Shannon Law states that the 
channel capacity, C is given by 



&=B\og*{l + {S/N)) b/s 



(17) 



For example, a 10kHz channel operating at a SNR of 15dB has a theoretical maximum information rate of 
10000 log 2 (3 1.623) = 498286A. 

METHODS 

Figure 1 the block design, it shows a model that has three main functions. It generates the position, velocity, and 
acceleration in polar (range-bearing) coordinates; it adds measurement noise to simulate inaccurate readings by 
the sensor; and it uses a Kalman filter to estimate position and velocity from the noisy measurements. 
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Figure l(a and b): Block diagram of the study design {Kalman filter tracking of the ComSat} 

The model is designed with the function kalman, having the following input data below 

1 . Kalman Filter 

1 . Range-bearing to XY 

2. XY Acceleration Model 



Table 1 : Showing Simulation Parameter 



Simulation Parameter 


Value 


Solver 


FixedStepDiscrete 


RelTol 


le-3 


AbsTol 


le-6 


Refine 


1 


MaxStep 


0.1 


MaxOrder 


5 


ZeroCross 


On 



Table 2: Band-Limited White Noise. Block Properties 



Name 


Cov 


Ts 


Seed 


Measurement Noise 


[2 2] 


0.1 


[52341 62341] 


Random satelitte Motion 


[0.001 .001] 


0.1 


[58413 74925] 





Table 


3: Gain Block Properties 












Name 


Gain 


Multiplication 


Param 
Min 


Param 
Max 


Param Data Type 
Str 


Out 
Min 


Out 
Max 


Out Data 
Type Str 


Meas. 
Noise 
Intensity 


[300 
0.01] 


Element- 
wise(K.*u) 


[] 


[] 


Inherit: Same as 
input 


[] 


[] 


Inherit: Same 
as input 



Table 4: Showing Block Types and Block Names 



BlockType 


Count 


Block Names 


Subsystem 


2 


Kalman Filter, XY Acceleration Model 


Outport 


2 


Est. Position [x, xdot, v, vdotl, Residuals 


Discretelntegrator 


2 


XY Position, XY Velocitv 


Band-Limited White Noise, 
(m) 


2 


Measurement Noise, Random satelitte Motion 


xy2rtheta (m) 


1 


XY to Range-Bearing 


Sum 


1 


Sum 


Gain 


1 


Meas. Noise Intensity 
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Table 5: Model Variables Showing Initial Value 



Variable Name 


Parent Blocks 


Calling string 


Value 


Speed 


XY Position 
XY Velocity 


[0,Speed] 
[0,Speed] 


400 



First we specify the satellite model with the process noise. 

Consider the discrete linear system with additive Gaussian noise w k on the input u k data. 



l *+l 



= Jx k + Btt k +w k 



ft = Cx k 



(18) 
(19) 



We want to use the available measurements y to estimate the state of the system x. We know how the system 
behaves according to the state equation, and we have measurements of the position, so how can we determine 
the best estimate of the state xl We want an estimator that gives an accurate estimate of the true state even 
though we cannot directly measure it. For this system w is the process noise and z is the measurement noise. We 
have to assume that the average value of w is zero and the average value of z is zero. We have to further assume 
that no correlation exists between w and z. having in mind the properties of a white noise. That is, at any time k, 
w h and Zk are independent random signal variables. Then the noise co variance matrices S w and S z are defined as 

Process noise co variance 



5. = £ KO 



Measurement noise covariance 



s, = £(-vf) 



(20) 



(21) 



where w T and z T indicate the transpose of the w and z random noise vectors, and E(-) means the expected value. 
There are many alternative but equivalent ways to express the Kalman filter equations. One of the formulations 
is given as follows: 



. ] 



** + i = (^t + But) + Kk{yt+i - ex*) 

P M = AP k A T + \. - APi&Sl'CPtA 7 



(22) 
(23) 
(24) 



The K matrix is called the Kalman gain, and the P matrix is called the estimation error covariance. The aim is to 
estimates the output y k given the inputs u k and the noisy output measurements 

>\- = C* t + h 

where zr is some Gaussian white noise. 

The first task during the measurement update is to compute the Kalman gain, The next step is to actually 
measure the process to obtain , and then to generate an a posteriori state estimate by incorporating the 
measurement. The final step is to obtain an a posteriori error covariance estimate. 

After each time and measurement update pair, the process is repeated with the previous a posteriori estimates 
used to project or predict the new a priori estimates. This recursive nature is one of the very appealing features 
of the Kalman filter — it makes practical implementations much more feasible than an implementation of a 
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Wiener filter which is designed to operate on all of the data directly for each estimate. The Kalman filter instead 
recursively conditions the current estimate on all of the past measurements. (Welch and Bishop 2006) 
The equations of the steady-state Kalman filter for this problem are given as follows. 

Measurement Update 

x (k) = x (k-l) + Pk(yk - C x (k-1) 



Time Update 



In these equations: 

• X(k-i) is the estimate of x (k) given past measurements up to y k _] 

• X(k) is the updated estimate based on the last measurement 

Given the current estimate, x (k) the time update predicts the state value at the next sample k+1 (one-step-ahead 
predictor). The measurement update then adjusts this prediction based on the new measurement y k+1 . P is the 
correction term is a function of the innovation, that is, the discrepancy.between the measured and predicted 
values of y k+1 . The innovation gain K k (Kalman gain) is chosen to minimize the steady-state covariance of the 
estimation error given the noise co variances. This filter generates an optimal estimate of y e . 

Note that the filter state is x (k . 1) 

Example: Given that 
A = 

3.1270 -0.9430 339 

1.0000 

1.0000 

B = 

-0.7660 

0.8190 

0.7330 

C = 

1 

D= 

Now we can construct a state-space model of the block diagram with the functions parallel and feedback. We 

build a complete mathematical model with u,w and as inputs y and y v (measurements) as outputs. 

a = A; 

b = [BB0*B]; 

c = [C;C]; 

d=[0 0;0 1]; 

P = ss(a,b,c,d,-l,'inputname',{'u' V V'},... 

'outputname',{y 'y v '}); 

To simulate the filter behavior we generate a sinusoidal input u and process and measurement noise vectors w 

and v. 

t = [0:100] time is from to lOOsecs, u = sin(t/15), Sw = 3.2, Sz = 1, n = length(t) =101 

randn('seed',0), w = sqrt(Sw)*randn(n,l), v = sqrt(Sz)*randn(n,l) 

we simulate to get 

[out,x] = lsim(SimModel,[w,v,u]); 

y = out(:,l); true response, ye = out(:, 2); filtered response, yv = y + v; measured response 
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The Kalman Gain K = 

1.0e+005 * 

4.5195 0.0000 0.0065 

0.0000 0.0000 0.0000 

0.0065 0.0000 0.0000 



Kalman filter response 



10 20 30 40 50 60 

No. of samples 



70 80 90 100 




10 20 30 40 50 60 70 80 90 100 
No. of samples 



Figure2: Kalman Filter Response (output and error) 

Measured Error is 

MeasErr = y-yv; 

Measured Error Covarianceis given by 

MeasErrCov = sum(MeasErr.*MeasErr)/length(MeasErr); 

Estimated Error is 

EstErr = y-ye; 

Estimated Error Covariance is given by 

EstErrCov = sum(EstErr.*EstErr)/length(EstErr); 

Measured Error Covariance is =0. 1657 

Estimated Error Covariance is =0.1 657 

Note since the measured error covariance and the estimated error covariance have same value we are sure of 

determining an unbiased estimate of the position of the Communication satellite. 

Simulation of Position of the Communication Satellite 

Initial Velocity Mismatch 

The Kalman Filter block works best when it has an accurate estimate of the satellite position and velocity, but if 

given time it can compensate for a bad initial estimate. 

To see this, 

Given Initial condition for estimated state parameter in the Kalman Filter as follows, 
Initial Satellite E-W estimate Position, with initial velocity of 400 in the Y direction, 
time t = [0 to 100] seconds at intervals of 10s, we have the result in figure 3 below: 
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Figure 3 Initial Satellite E-W estimate Position, with initial velocity of 400 in the Y direction 

RESULTS AND ANALYSIS 

Estimation of the position and velocity is performed by the 'Kalman Filter' subsystem. This subsystem samples 
the noisy measurements, converts them to rectangular coordinates, and sends them as input to the Signal 
Processing block set, Figure lb 
The Kalman Filter block produces two outputs in this application. 

• An estimate of the actual position. This output is converted back to polar coordinates so it can be 
compared with the measurement to produce a residual,or error( the difference between the estimate 
and the measurement). The Kalman Filter block smoothes the measured position data to produce its 
estimate of the actual position. 

• The second output from the Kalman Filter block is the estimate of the state of the satellite, In this case, 
the state is comprised of four numbers that represent position and velocity in the X and Y coordinates. 

Changing the estimate from 400 to 100 and running the model again, we have 
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Figure 4: Showing Velocity in Y direction=100 

It is observed that the range residual is much greater and the 'E-W Position' estimate is inaccurate at first. 
Gradually, the residual becomes smaller and the position becomes more accurate as more measurements are 
gathered. 

Simulation by increasing the Measurement Noise Estimate 

In the present model, the noise added to the range estimate is rather small compared to the ultimate range. The 
maximum magnitude of the noise is 300 ft, compared to a maximum range of 40,000 ft. Increasing the 
magnitude of the range noise to a larger value, for example, 5 times this amount or 1500 ft. by changing the first 
component of the Gain parameter in the 'Meas. Noise Intensity' Gain block. 
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Figure 5: Poor Prediction 

We observe that the blue lines representing the estimated positions have moved farther from the red lines 
representing the actual positions, and the curves have become much more 'bumpy' and 'jagged'. We can partially 
compensate for the inaccuracy by giving the Kalman Filter block a better estimate of the measurement noise. 
Setting the Measurement noise covariance parameter of the Kalman Filter block to 1500 and running the model 
again. 
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Figure 6 Showing Unbiased Estimate 

ANALYSIS AND DISCUSSION 

• From Figure 2 the covariance plot, one can see that the output covariance did indeed reach a steady 
state in about five samples., the Co variances arrived at the same value, being the same, that means that 
estimated output of our kalman filter is unbiased, thereby achieving the aim of the study. 

• From figure 5, It is noted that when the measurement noise estimate is better, the E-W and N-S 
position estimate curves become smoother. This is the expected behavior. 

CONCLUSION AND RECOMMENDATION 

Conclusion: For communication satellites launched in the orbit, the intricacy of predicting its position in the 
orbit is high due to the effect of gravity waves, drag resistance, solar arrays etc, which give rise to erroneous 
range detection of the space craft. The kalman filter has profer a solution to the mentioned problems, it is 
important to state here that it not only works well in practice, but it is theoretically attractive because it can be 
shown that of all possible filters, it is the one that minimizes the variance of the estimation error. 

RECOMMENDATION 

It is recommended that the Kalman filter is being used to filter noise in order to track or predict linear systems. 
Many physical processes, can be approximated as linear systems and can be predicted e.g. vehicle driving along 
a road, a motor shaft driven by winding currents, or a sinusoidal radio-frequency carrier signal. 
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